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Abstract. We calculate the exact stationary distribution of the one-dimensional 
zero-range process with open boundaries for arbitrary bulk and boundary hopping 
rates. When such a distribution exists, the steady state has no correlations between 
sites and is uniquely characterized by a space-dependent fugacity which is a function 
of the boundary rates and the hopping asymmetry. For strong boundary drive the 
system has no stationary distribution. In systems which on a ring geometry allow 
for a condensation transition, a condensate develops at one or both boundary sites. 
On all other sites the particle distribution approaches a product measure with the 
finite critical density pc- In systems which do not support condensation on a ring, 
strong boundary drive leads to a condensate at the boundary. However, in this case 
the local particle density in the interior exhibits a complex algebraic growth in time. 
We calculate the bulk and boundary growth exponents as a function of the system 
parameters. 
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1. Introduction 

The zero-range process (ZRP) has originally been introduced in 1970 by Spitzer [1] as a 
model system for interacting random walks, where particles on a lattice hop randomly to 
other neighboring sites. The hopping rates Wn depend only on the number of particles 
n at the departure site. This model has received renewed attention because of the 
occurrence of a condensation transition [2-6] analogous to Bose-Einstein condensation 
and because of its close relationship with exclusion processes [7]. Condensation 
phenomena are well-known in colloidal and granular systems (see [8] for a recent study 
making a connection with the zero-range process), but appear also in other contexts, 
such as socio-economics [9] and biological systems [10] as well as in traffic flow [11] 
and network theory [12, 13]. In the mapping of the ZRP to exclusion processes in 
one space dimension, condensation corresponds to phase separation. The ZRP has 
served for deriving a quantitative criterion for the existence of non-equilibrium phase 
separation [14] in the otherwise not yet well-understood driven diffusive systems with 
two conservation laws. For the occurrence of condensation the dimensionality of the 
ZRP does not play a role and we shall consider only the one-dimensional (Id) case. 

Most studies of the ZRP focus on periodic boundary conditions or the infinite 
system. Under certain conditions on the rates (sec below) the grand-canonical 
stationary distribution is a product measure, i.e. there are no correlations between 
different sites [15]. In addition, an exact coarse-grained description of the dynamics 
is possible in this case in terms of a hydrodynamic equation for the particle density 
p{x,t) [16,17]. 

In zero-range processes for which the hopping rates w„ admit condensation, one 
finds that above a critical density pc a finite fraction of all particles in the system 
accumulate at a randomly selected site, whereas all other sites have an average density 
Pc [5,18, 19]. The large scale dynamics of condensation has been studied in terms of a 
coarsening process [19,20]. The steady-state and the dynamics of open systems which 
may admit condensation has not been addressed so far. 

In the present work we consider a ZRP on an open chain with arbitrary hopping 
rates and boundary parameters. Particles are added and removed through the 
boundaries. In the interior of the system hopping may either be symmetric or biased 
in one direction. We calculate the exact steady-state distribution, when it exists, and 
find it to be a product measure. In order to study condensation phenomena in the open 
system we analyze in detail a particular but generic ZRP which admits condensation. 
In this model the hopping rates for large n take the form w„ = 1 + 6/n. On a periodic 
lattice it is known that a condensation takes place at high densities when b > 2. We 
find that in an open system and for a weak boundary drive the model evolves to a non- 
critical steady-state. On the other hand, if the boundary drive is sufficiently strong, 
the system may develop a condensate on one or both of its boundary sites even for 
b < 2. The number of particles in the condensate grows linearly in time due to the 
strong boundary drive. In this case the interior of the system may either (a) reach a 
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Figure 1. Graphic representation of the ZRP model with open boundaries. 



sub-critical steady-state, (b) reach a critical steady-state, or (c) it may evolve such that 
the local particle density exhibits a complex algebraic growth in time. Which behavior 
is actually realized depends on the boundary rates, on the parameter b, and on the 
asymmetry in the hopping rates. 

The paper is organized as follows. In Sec. |2l we define the ZRP with open 
boundaries, and give some examples of possible hopping rates Wn- The exact steady- 
state distribution for general boundary parameters and rates Wn is derived in Sec. EJ 
In Sec. m we consider the behavior of finite systems for totally asymmetric, partially 
asymmetric and symmetric hopping rates. The long-time temporal behavior of local 
densities is obtained, and supporting numerical simulations are presented. The long- 
time behavior of bulk densities in an infinite system is exactly obtained using a 
hydrodynamical approach in Sec. |S| Finally, we present our summary and conclusions 
in Sec. IHI 

2. Zero-range processes with open boundaries 

The ZRP on an open Id lattice with L sites is defined as follows. Each site k may be 
occupied by an arbitrary number n of particles. In the bulk a particle at site k (say, 
the topmost of n particles) hops randomly (with exponentially distributed waiting time) 
with rate pWn to the right and with rate qWn to the left. Without loss of generality, we 
take throughout this paper p > q, so that particles in the bulk are driven to the right. 
At the boundaries these rules are modified. At site 1 a particle is injected with rate a, 
hops to the right with rate pWn, and is removed with rate •ywn- At site L a particle is 
injected with rate 6, hops to the left with rate qWn, and is removed with rate (see 
figure HI) . 

Examples for such processes include the trivial case of non-interacting particles 
{wn = n) and the pure chipping process (w„ = 1) [21]. In the mapping to exclusion 
processes where the particle occupation number becomes the interparticle distance, the 
chipping model defined on a ring maps onto the simple exclusion process which is 
integrable and can be solved by the Bethe ansatz, combinatorial methods, recursion 
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relations and matrix product methods [22,23]. A largely unexplored but interesting 
integrable model with rates w„ = (1 — — q) that interpolates between these 

two cases has been found in [24]. The truncated chipping process with w„ = 1 for 
1 < n < K and Wn = oo for n > K maps onto the drop-push model [25] which is 
also integrable [26] . For wi = w and Wn = ^ for n > 2 one obtains the version of the 
non-integrable KLS model [27] that has been introduced in [28] as a toy model for traffic 
flow with a nonsymmetric current density relation. A parallel updating version of this 
model was found to correspond to a broader class of traffic models [29] as well as an 
integrable family of ZRP's [30] which includes the asymmetric avalanche process [31]. 

In a similar mapping to exclusion processes (but in a continuum limit) the family 
of models with hopping rate 

(1) 

describes the diffusion of interacting rods on the real line [32]. The condensation 
model [4, 5, 14, 18-20] with 

b 

Wr,= l + - 2 

n 

is a generic model that exhibits the condensation phenomenon described in the 
introduction for b > 2. Both models are non-integrable and hence alternative tools 
must be employed for deriving information about their dynamical behavior. We remind 
the reader that this model is generic in the sense that it represents the complete family 
of models with rates of the form Wn — 1 + b/n + 0{rr^) where s > 1. Other ZRP's 
which exhibit condensation are defined by the rates Wn — l-\-b/rf with < cr < 1 and 
6 > [14] or by rates approaching zero [6]. 



Wr. 



1 + ^ 

n 



3. Stationeiry distribution 

A state of the model at time t may be defined through a probability measure P„ on 
the set of all configurations n — {ni,n2, ■ ■ ■ jTIl), e N. Here is the number of 
particles on site k. To calculate the stationary distribution it is convenient to represent 
the generator H of this process in terms of the quantum Hamiltonian formalism [23] 
where one assigns a basis vector \n) of the vector space to each configuration 

and the probability vector is defined by \P) = X^n-^^N)- normalized such that 

{s\P) = 1 where (s| = {n\n') = 5n,n'- The time evolution described above 

is given by the master equation 

||P(t)) = -H\P{t)) (3) 

through the "quantum Hamiltonian" H. This operator has off-diagonal matrix elements 
Hn^n' which are the hopping rates between configurations n, n' and complementary 
diagonal elements to preserve conservation of probability. 
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Since we have only nearest neighbor exchange processes the Hamiltonian in (0) can 
be written as 

L-l 

H = hi + hL + ^hk,k+i, (4) 
fc=i 

where hk^k+i acts nontrivially only on sites k and k + 1 (corresponding to hopping) 
while hi, hi generates the boundary processes specified above. For the ZRP we define 
the infinite-dimensional particle creation and annihilation matrices 
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as well as the diagonal matrix d with elements dij 
have 



Wi6ij. With these matrices we 



hk. 



k+l 



k+l 



d, 



k+l) 



and 



hi = a{a^ — 1) + 7(^1 — di), —h^ = S{a^ — 1) + ~ ^l)- 



(6) 

(7) 



The "ground state" of H has eigenvalue 0. The corresponding right eigenvector is the 
stationary distribution which we wish to calculate. 

Guided by the grand-canonical stationary distribution of the periodic system we 
consider the grand-canonical single-site particle distribution where the probability to 
find n particles on site k is given by 



PHrik 



n] 



n 



(8) 



Here the empty product n = is defined to be equal to 1 and Zk is the local analogue 
of the grand-canonical partition function 



Zk^Z{zk) = Y.zl\{ 



(9) 



n=0 



i=l 



The corresponding probability vector \P^) with the components P*{nk = n) satisfies 



^-k'd\P*k, 



a'\Pt 



^i^Pk) 



(10) 



The proof of this property is by straightforward calculation. 

As an ansatz for calculating the stationary distribution we take the L-site product 
measure with the one-site marginals (jHl) which is given by the tensor product 



P*) = \P* 



\Po 



\PI) 



(11) 
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and according to (jTU)) satisfies 



H\P*) 



L-l 
k=l 

+ {a- iz,){z^H, - 1) + (5 - Pzl){ziHl - 1)] I P* ). (12) 

The uncorrelated particle distribution (lllj) is stationary if and only if all terms on the 
right hand side of this equation cancel. It is not difficult to see that this is satisfied for 
the following stationarity conditions on the fugacities z^ 

pzk - qzk+i = a--fzi = (3zl -6 = c . (13) 

The quantity c is the stationary current. 

This recursion relation has the unique solution 

[(« + 6){p-q)-af3 + 75] (l)'"' - 7<5 + (f)^"' 
lip-q-P)+P{p-q + 7) (f ) 

and the current is given by 

-75 + a/5 

c={p-q) ^ — 7-T- (15) 



l{p - q - (3) + (3{p - q + 7) 



^9, 

Thus, given the rates Wn the stationary distribution is unique and completely specified 
by the hopping asymmetry and the boundary parameters. The stationary density profile 
follows from the fugacity profile through the standard relation 

d 

Pk = ZkT:— In Zk , (16) 

OZk 

where Zk = Z{zk) is determined by the bulk hopping rules, as given in (0). 

We remark that up to corrections exponentially small in the system size L, the bulk 
fugacity is a constant 

ZeS = — , (17) 

p - g + 7 

that depends only on the left boundary rates. This is in agreement with the observation 
[17] for special boundary parameters and can be explained in terms of the more general 
theory of boundary- induced phase transitions [33] . 

Some special cases of (fTHl deserve mentioning. 
(1) Setting the boundary extraction rates equal to the corresponding bulk jump rates, 
i.e., P = p, 'J = q, and defining reservoir fugacities Zr^i by a = pzi, 6 = qzr, the 
expression (fT^ reduces to 

~ / \L-l 

1_ fP 
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This dynamics has a natural interpretation as couphng the system to boundary reservoirs 
with fugacities Zr, zi respectively. This case has been considered in [17] and earlier 
in [15,34] for p = q. The general solution ()14|1 for arbitrary boundary rates appears to 
be a new result for the ZRP. 

(2) For 

^S = a(3(^] (19) 

the current vanishes and the system is in thermal equilibrium. For the symmetric case 
p = q this implies constant fugacities Zk- 

(3) For symmetric hopping p = q = 1 the fugacity profile is generally linear 

_ a + 6 + a(3{L - 1) - {a(3 - -f6){k - 1) 

P + l + Pl{L-l) ' ^''^ 

with a current 

a/3 — 7(5 



^ /3 + 7 + /37(L-l) • ^^^^ 
Notice that the linear fugacity profile does not imply a linear density profile except in 
the very special case of non-interacting particles where Zk = Pk- 
(4) In the totally asymmetric case, g = 0, we find 

Zk = = zioik^ L , (22) 

The current is given hy c = pz. 



In the considerations above we have tacitly assumed that the local partition function 
Zk exists for all k. Since for suitable choices of boundary parameters the local fugacities 
at the boundary may take arbitrary values this assumption implies an infinite radius of 
convergence of Zk = Z{zk) which is not the case in the models described above. This 
raises the question of the long-time behaviour of the ZRP for strong boundary drive, 
i.e., rates which drive the boundary fugacities out the radius of convergence of Z. This 
is studied in detail in the following sections for the condensation model The rod 
model dH) can be regarded as a generic ZRP without bulk condensation transition, but 
finite radius of convergence for Z. Where appropriate we compare its behaviour with 
that of the condensation model. 



4. Condensation model — Steady state and dynamics near the boundary 

In this section we consider in some detail the long-time behavior of the condensation 
model ((21). For this model the local fugacity Zk has to satisfy < 1 in the steady state. 
On a ring geometry, the model exhibits a condensation for 6 > 2 at high density. We 
first analyze the totally asymmetric case, where particles can only hop to the right. The 
cases of partial asymmetry and symmetric hopping are then treated. 
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4.1. Totally asymmetric hopping 

For the totally asymmetric case we take g = 7 = 5 = 0. For the normalization of time 
we set p = 1. The exact steady-state solution (fT^ yields 

Zk = a = z ioT k L , (24) 

ZL = ^, (25) 

and the current is given hj c = z. Since for this model z has to satisfy z < 1, the steady 
state (j211l2Sl) is valid only for a < 1 and P > a. In this case the single-site steady-state 
distribution is given, for large n, by P*{nk = n) ~ z"'/n^. 

We now proceed to discuss the dynamical behavior of the model in the case that 
stationary state does not exist. For a > 1 the following picture emerges. 

Site 1: 



On site 1 particles are deposited randomly with rate a > 1 and are removed by hopping 
to site 2 with rate 1 + b/ni. Hence the occupation number performs a simple biased 
random walk on the set ni of positive integers with drift a — 1 — b/ni which is positive 
for any ni > b/ (a — 1). Such a random walk is non-recurrent and reaches asymptotically 
the mean velocity v = a — 1. Hence the mean particle number Ni{t) = ( ni(t) ) on site 1 
grows asymptotically linearly 

iVi(t) ~ (a - l)t. (26) 



Boundary sites k > 1: 



We extend the random walk picture (which is strictly valid for site 1) to site 2. On site 2 
particles are injected (by hopping from site 1) with rate 1 + fo/rii and are removed with 
rate 1 + b/n2. Since ui increases in time on average the input rate approaches 1 and 
the occupation number at site 2 performs a biased random walk with hopping rate 1 to 
the right and rate 1 + 6/n2 to the left. Whether this random walk is recurrent depends 
on b. The asymptotic behavior has been analyzed in [35]. We merely quote the result: 



N2{t) ~ < 



^1/2 

t^/yint 

\nt 

P* 



b< 1 
b= 1 
l<b<2 
b = 2 
b>2 



(27) 



The constant 



P 



6-2 



f28) 



is the critical density of the condensation model [4]. It is approached with a power 
law correction By applying this random walk picture to further neighboring 

sites, and assuming scaling, it has been shown that neighboring boundary sites behave 
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Figure 2. Temporal evolution of local densities, as obtained from numerical 
simulations of the totally asymmetric condensation model, with b = 3/2, and a — 2. 
The solid line corresponds to the expected growth law t^^^, and the dotted line 
corresponds to linear growth. 



asymptotically in the same fashion [35]. Similar analysis shows that the square-root 
increase of the particle density occurs also for the model ((T)) for all values of its interaction 
parameter b'. 

To test the validity of the random walk picture to sites beyond = 2, we carried 
out numerical simulations of a totally asymmetric model with L = 5. We first note that 
since the hopping is totally asymmetric the time evolution of the system up to site k 
is independent of what happens at sites to the right of k. In particular, the dynamics 
on all sites k < L is independent of /3. Hence, in order to study boundary layers it is 
sufficient to simulate very small systems of only a few sites. In figure El we present the 
long-time behavior of the occupation of sites k = 1 — 4 for a = 2 and b = 3/2. It is 
readily seen that while site 1 grows linearly in time, sites 2 — 4 grow with the expected 
power law t^^'^. 

Bulk sites A; > 1: 

The picture of the simple random walk becomes increasingly inaccurate as one enters into 
the bulk of the system, since injection events onto a site become increasingly correlated 
in time. This violates the random walk assumption and makes the previous analysis 
invalid in this case. The temporal behavior of bulk sites k ^ 1 can be treated exactly by 
using a hydrodynamic approach, which yields the behavior of bulk sites in the long-time 
limit. This analysis, carried out in Sec. shows a different dynamical behavior in the 
bulk. Notice, however, that for any finite system all "bulk sites" have finite distance 
from the boundary and hence behave asymptotically like the boundary sites. 

Site L: 

Since the motion of particles on site k is independent of the motion on sites to their 
right we expect the bulk result to be asymptotically valid on all sites up to site L — 1, 
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Figure 3. Temporal evolution of local density at the rightmost site, as obtained from 
numerical simulations of the totally asymmetric condensation model, with b = 3/2, 
a = 2 and (3 = 1/2, 1, 2. The solid line corresponds to the expected growth law t^^^, 
and the dotted line corresponds to linear growth. 



i.e., there is no right boundary layer with yet another set of growth exponents. On site 
L the following picture emerges. There is an asymptotic incoming flux c = 1 and an 
exit rate (3{1 + b/ni)- For f3 = 1 the boundary site behaves like a bulk site and we 
obtain the bulk growth exponent. For (3 < 1 the outgoing flux does not compensate the 
incoming flux which yields asymptotically linear growth piit) = (1 — l3)t. For /3 > 1 a 
finite stationary chemical potential = 1//3 is approached. 

In figure El we present simulation results for the long-time dynamics of the 
occupation of site k = L = 5. Depending on the value of j3 the occupation number 
either grows linearly (/3 < 1), grows with the same power law as the bulk {j3 = 1), or 
approaches a finite density (/3 > 1), as expected from the above discussion. 

To complete the discussion of the totally asymmetric case we remark that for a 
subcritical left boundary a < 1 but supercritical P < a the bulk of the system becomes 
stationary with fugacity Zbuik = On site L a condensate develops with linearly 
increasing particle density pL(t) ~ (a — (3)1. 

4-2. Partially asymmetric hopping 

We now analyze the dynamical behavior of the partially asymmetric model in the case 
where no stationary state exists. We start by considering the case where the rates at 
both boundaries are such that the occupation of the two boundary sites increase linearly 
with time. This takes place for a — •y < p — q < 13 — 5. Here sites k = 1 and k = L act as 
reservoirs for the rest of the system. The effective rates at which the reservoirs exchange 
particles with the system are aes = PeS = P and 7efr = ^efr = Q- The fugacity ((Hj) at 
sites k = 2, ... ,L — 1 thus becomes = 1. The asymptotic temporal behavior ()27|) 
holds also in this case. 

The picture changes qualitatively if only one boundary fugacity does not exist. 
Suppose first that this happens at site 1. In this case the density on site 1 will increase 
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Figure 4. The fugacity profile H14|) for the partiaUy asymmetric case, witli b = 3/2, 
p — 3/A and (7 = 1/4. Solid line corresponds to a left condensate (a = 3/4,/3 = 1,7 = 
1/4, S = 1/5) and dashed line corresponds to a right condensate {a = 1/2, f3 — 3/4, 7 = 
1/4,(5= 1/4). 

linearly with time, as in the previous totally asymmetric case. This happens when 
a — 7 > p — q and (3 — 5 > p — q. As before, site 1 acts as a reservoir for the rest 
of the system, with effective rates a^s = p, 7efT = q- Sites k = 2, . . . , L — 1 are thus 
stationary, with the fugacity given by (fT^ . The fugacity at sites away from the right 
boundary approaches 1, with a deviation exponentially small in the system size (see 
figure 131). Hence, starting from an empty initial state one expects algebraic growth of 
the local density until, after a long crossover time, stationarity is reached. 

If, on the other hand, the right boundary rates drive site L out of equilibrium, then 
the density on site L increases linearly in time. This happens when a — 'y < p — q and 
(3 — 5 < p — q. Site L acts effectively as a boundary reservoir with /Seg = p and 5eff = q- 
As in the preceding case the system becomes stationary, but with a finite (independent 
of system size) deviation of the bulk fugacity from the critical value z = 1. Both in the 
bulk and at the left boundary the density approaches the finite value dictated by the 
left boundary fugacity (see figure Ej). 

Finally, for (3 — 6<p — q<a — both boundary sites have finite fugacity, and the 
system reaches a non-critical steady state, as given by p4p . 

4-3. Symmetric hopping 

The analysis of symmetric hopping {p = q = 1) follows very closely that of the partially 
asymmetric case, except that here the fugacity profile is linear rather than exponential. 
In particular, if only one boundary is driven out of equilibrium, a condensate appears 
at that boundary with a particle density which increases linearly with time. This site, 
say site 1, acts as a reservoir with Ues = Jes = 1- The fugacity profile decreases linearly 
from 1 — 0(1/L) at the left boundary to 6/ P + 0(1/ L) at the right boundary, as given 
by (j2ni)- 

When both boundary sites are supercritical they act as reservoirs with effective 



left condensate 



right condensate 
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rates a^s = f^es = 7efr = (^efr = 1, yielding Zk = I throughout the system. This is 
also the case relevant for studying the fluid-mediated interaction of probe particles in 
two-species systems [35]. We expect similar boundary growth laws as in the asymmetric 
case. This is because the incoming flux from the direction of the bulk is compensated 
by the outgoing flux at the boundary and one has an effective random picture with the 
same rates as above. However, corrections to scaling are expected to be larger as current 
correlations are now stronger due to hopping contributions from the bulk. Moreover, 
as opposed to the asymmetric case, in a semi-infinite system both boundaries have the 
same growth exponents. 



5. Hydrodynamics — Exact analysis of dynamics in the bulk 

In order to analyze the time evolution of sites far away from the boundaries we consider 
the hydrodynamic limit of the ZRP model. The coarse-grained time evolution of the 
density profile starting from a non-stationary initial profile can be determined, by 
adapting standard arguments [16,36], from the continuum limit of the lattice continuity 
equation. This equation reads 
d 

—Pk = Ck-i - Ck , (29) 
with the local current 

Ck = pzk - qzk+i , (30) 

and z expressed in terms of p. Together with an appropriate choice of constant boundary 
fugacities the solution is uniquely determined [17]. 

For the driven system one obtains under Eulerian scaling (lattice constant a ^ 0, 
t ^ t/a, system length fixed) 

d d 

— p(x, t) = -{p - q)—z{x, t) + a{p + q)/2—z{x, t) (31) 

where the infinitesimal viscosity term serves as regularization and selects the physical 
solution of the otherwise ill-posed initial value problem with fixed boundary fugacities. It 
takes care of a proper description of discontinuities which may arise in the form of shocks 
or a boundary discontinuity. Indeed, in the large-time limit the density approaches a 
constant given by the left boundary fugacity, with a jump discontinuity at the right 
boundary [17]. This is in agreement with the exact result derived in the previous 
subsection. We stress that (p?T|) provides an exact description of the density evolution 
under Eulerian scaling. It is not a continuum approximation involving a mean field or 
other assumption. The simple form of (j!?T|l originates in local stationarity. The absence 
of noise is due to Eulerian scaling, i.e., the effects of noise appear on finer scales. For 
a recent rigorous discussion of the hydrodynamic limit of stochastic particle systems, 
see [16,37]. 

Consider a semi-infinite system which is initially empty and has constant density 
at the left boundary. This boundary condition induces a rarefaction wave entering the 
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bulk which can be constructed using the method of characteristics. The speed Vq of the 
wave front is given by the zero-density characteristic of ()31|) which is the average speed 
■^0 = {p — of a single particle. In light of the results of the previous sections we are 
particularly interested in the case where the left boundary fugacity is equal to 1. For 
6 < 2 in the condensation model [and for any h' in the rod model (^] this corresponds to 
an infinite boundary density. On the other hand, for h > 2 the corresponding boundary 
density is pc (EH)- Therefore we are searching for a scaling solution in terms of the 
scahng variable u = x/iyot) such that p{u) = for m > 1 while p = pc or p = oo at 
the left boundary. On physical grounds the solution has to be continuous as no shock 
discontinuities can develop for the initial state (empty lattice) under consideration. 
Under these conditions (jHTj) can be integrated straightforwardly by setting a = and 
one finds the implicit representation 
dz 

— = uwi . (32) 

dp 

This is the solution within the interval Vit < x < VqI. Here 

dz 

Vz = {p-q)-r- (33) 
dp 

is the collective velocity of the lattice gas which plays the role of the speed of the 
characteristics for the hydrodynamic equation 1)3111 . Outside this interval one has z = 
for X > VqI and z = 1 for < x < Vit. We now analyze this solution in terms of p. 

For the model (0) one has Z = 1/(1 — z)^'^^ which yields the fugacity-density 
relation z = p/ {b' + 1 + p). Hence 



piu) = {b' + l){J--l]. (34) 




with Vq = {p — q)/{b' + l) and f i = for all b'. At any fixed x the bulk density increases 
algebraically 

Pbuikit) ~ t^/^ (35) 

for any 6' > 0. 

For the condensation model the local partition is the hypergeometric function 

Z = 2Fiil,l-A + b;z) (36) 

which does not admit an explicit representation of 2; as a function of p. However, as the 
fugacity is an increasing function in time approaching 1 we may analyze its asymptotic 
behavior by expanding the hypergeometric function around z = 1 [19]. 

b<2: 

Here the left boundary fugacity z = 1 corresponds to infinite left boundary density. 
Moreover vi = 0, therefore the solution ()32p describes the density profile in the interval 
< u < 1. This enables us to consider fixed x and study the long time limit. In order 
to analyze the small u behavior (i.e. z close to 1 and p large) one has to distinguish two 
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domains [19]. For 6 < 1 one has p{z) oc z/{l — z) for large p while for 1 < 6 < 2 one 
finds p{z) (X z/{l — z)"^"^. As a function of u we make the ansatz p oc for the large 
t (i.e. small u) asymptotics. Inserting this into the differential equation (|!?T|l yields a 
consistent solution with 

Phuik{t) ~ e (37) 

only for 

We conclude that the bulk density increases algebraically with the universal diffusive 
exponent 1/2 for 6 < 1 and with a non- universal 6-dependent exponent in the range 
1 < 6 < 2 of the condensation model. At h = 1,2 there are logarithmic corrections 
which we do not discuss further. 

For symmetric hopping one describes the density dynamics under diffusive scaling 
t — > t/a^ and obtains 

d d'^ 

which needs no further regularization. Repeating the previous analysis one finds the 
same bulk growth exponents as for the asymmetric hopping model. 



b>2: 

In the condensation regime the boundary density corresponding to z = 1 is 
Pc (EHl)- Again two different regimes are found from the asymptotic analysis of the 
hypergeometric function. For 6 < 3 one has f i = and repeating the same analysis 
as for 1 < 6 < 2 shows that at fixed x the density approaches pc with a power law 
correction with the same exponent (2 — 6)/(3 — b) as before. This is analogous to the 
behavior in the boundary sites analysed in the previous section, but the exponent is 
different. For 6 > 3 one finds for the collective velocity [19] 

,(6-3)2(6-2)2 ^ , 

vi = {p-q y ^ >0- (40) 

Hence an analysis of the long-time behaviour at fixed x is not meaningful. A domain 
with constant p = Pc spreads into the system, with a front speed v = vi. This front is 
preceded by the rarefaction wave ()32p for Vit < x < v^t. This behaviour as a function 
of b is unexpected as usually changes in the rarefaction wave of this type are caused 
by changing the boundary density rather than an interaction parameter of the driven 
system. 

We stress that there are two questions that cannot be answered by the 
hydrodynamic analysis given above. First we apply Eulerian or diffusive scaling 
respectively. This leaves generally open what happens in any semi-infinite lattice system 
at finite lattice distance from the boundaries or in a finite system. Any deviation from 
the results given above which decays on lattice scale as one approaches the bulk cannot 
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be detected within the hydrodynamic description. Boundary layers, which have been 
analyzed in the previous section and found to have an interesting microscopic structure, 
would under scaling at best appear as a structureless boundary discontinuity. 

Secondly, the radius of convergence of Z is 1 and forcing the boundary fugacities 
zi or zl to be larger than 1 implies a breakdown of the assumption of local stationarity 
underlying the hydrodynamic description, at least in the vicinity of the boundaries. In 
particular, the hydrodynamic approach is not applicable for analyzing the condensation 
regime 6 > 2 if the boundary density exceeds the critical density of the bulk. However, 
this regime may be analyzed by the random-walk approach discussed in the previous 
section. 

6. Conclusions 

In this paper we studied the dynamical behavior of the zero-range process with open 
boundary conditions for arbitrary bulk and boundary rates. It is found that for a weak 
boundary drive the model reaches a steady state. The exact steady-state distribution 
is calculated and is shown to be a product measure characterized by site-dependent 
fugacities. In the case of strong drive the system does not reach a steady-state, and its 
evolution in time is calculated. To this end we considered the condensation model with 
an initially empty lattice and studied how the local density evolves in time. As long as 
the bulk dynamics does not permit condensation (6 < 2) the growth of the local density 
is algebraic in time with exponents that we determined using a random walk picture for 
the boundary region, and standard hydrodynamic description in the bulk. Prom this 
analysis we are led to the conclusion that in the condensation model with b > 2 (where 
condensation in a periodic system sets in above the critical bulk density pc — 1/(6 — 2)) 
only the boundary sites develop into a condensate, with a density increasing linearly 
in time. All bulk sites become stationary with finite local fugacities determined by the 
boundary rates. Somewhat surprisingly the driven and the symmetric model have the 
same bulk and boundary growth exponents. The boundary condensate appears also for 
b < 2 when no condensation transition exists in a periodic system. It is a general feature 
of zero-range processes in which Z has a finite radius of convergence. 

For 6 > 1 we observe a precursor to the condensation transition in the sense that 
the universal diffusive growth for & < 1 breaks down. Bulk and boundary growth 
exponents become different and both decrease with b, i.e. they become non-universal. 
It is interesting to note that 6 = 1 plays a special role also for the stationary state on a 
ring geometry. For b < 1 the stationary probability to find any given site empty vanishes 
as the density is increased to infinity [19], in agreement with intuition. However, for 
6 > 1 every site has a finite probability of being empty, even if the particle density is 
infinite. Applied to present scenario this implies the counterintuitive result that even at 
very late time, when the average particle density in an open system tends to infinity on 
each site, one still expects to find a finite fraction of empty sites at any given moment. 

In this context it is also instructive to study the mean first passage time (MFPT) 
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of the boundary particle density, i.e. the mean time r after which a particle number N 
has been reached for the first time at a given site, starting from an empty site. Using 
the exact general MFPT expression [38] for the effective random walk defined above one 
finds 



where for non-integer b the factorials are defined by the F-function. For large N this 
quantity has the asymptotic behavior 



Again there is a transition at 6 = 1, with diffusive behavior for b < 1 and sub-diffusive 
exploration of the state space for 6 > 1. In the bulk a simple random picture for 
the on-site density dynamics is not valid and a prediction for the bulk MFPT is not 
possible. For small finite system size one expects boundary behavior everywhere, but 
with increasing corrections to scaling as one moves away from the boundary. The precise 
nature of the crossover from the boundary growth exponent to the hydrodynamic bulk 
growth exponent remains an open problem. 
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